Space-time crystals of trapped ions 
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Great progresses have been made in exploring exciting physics of low dimensional materials in 
last few decades. Important examples include the discovering and synthesizing of fullerenes (zero- 
dimensional, 0D) [U, carbon nanotubes (ID) 0] and graphene (2D) A fundamental question 
is whether we can create materials with dimensions higher than that of conventional 3D crystals, 
for example, a 4D crystal that has periodic structures in both space and time. Here we propose 
a space-time crystal of trapped ions and a method to realize it experimentally by confining ions 
in a ring-shaped trapping potential with a static magnetic field. The ions spontaneously form a 
spatial ring crystal due to Coulomb repulsion. This ion crystal can rotate persistently at the lowest 
quantum energy state in magnetic fields with fractional fluxes. The persistent rotation of trapped 
ions produces the temporal order, leading to the formation of a space-time crystal. We show that 
these space-time crystals are robust for direct experimental observation. The proposed space-time 
crystals of trapped ions provide a new dimension for exploring many-body physics and emerging 
properties of matter. 



Symmetry breaking plays profound roles in many-body 
physics and particle physics [J] . The spontaneous break- 
ing of continuous spatial translation symmetry to discrete 
spatial translation symmetry leads to the formation of 
various crystals in our everyday life. Similarly, the spon- 
taneously breaking of time translational symmetry can 
lead to the formation of a time crystal 0, @] ■ Intuitively, 
if a spatially ordered system rotates persistently in the 
lowest energy state, the system will reproduce itself pe- 
riodically in time, forming a time crystal in analog of an 
ordinary crystal. Such a system looks like a perpetual 
motion machine and may seem implausible in the first 
glance. On the other hand, it has been known that a su- 
perconductor @,IH or even a normal metal ring [§- 11 1 can 
support persistent currents in its quantum ground state 
under right conditions. However, the rotating Cooper 
pairs or electrons in a metal are homogenous in space. 
Thus no spatial order is mapped into temporal order and 
it is not a time crystal Q. While it has been proved 
mathematically that time crystals can exist in principle 
@, 0] , h was not clear how to realize and observe time 
crystals experimentally. 

In this paper, we propose a method to create a space- 
time crystal (Fig. la), which is also a time crystal, with 
cold ions in a cylindrically symmetric trapping potential. 
Different from electrons in conventional materials, ions 
trapped in vacuum have strong Coulomb repulsion be- 
tween each other and have internal atomic states. The 



strong Coulomb repulsion between ions enables the spon- 
taneous breaking of the spatial translation symmetry, re- 
sulting in the formation of a spatial order that can be 
mapped into temporal order. The internal atomic states 
of ions can be utilized to cool the ions to the ground state 
as well as observing their persistent rotation directly by 
state-dependent fluorescence. 

Trapped ion Coulomb crystals have provided unique 
opportunities for studying many-body phase transitions 
12, 1 1 31 ] and quantum information science 0E3. The 
two most common types of ion traps are the Penning 
trap [3] and the Paul trap [lf| . The magnetron motion 
of an ion in the Penning trap is an orbit around the top 
of a potential hill that is not suitable for studying the 
behavior of the ion in its rotational ground state. Here we 
propose to use a combination of a ring-shaped trapping 
potential from a variation of the Paul trap, and a weak 
static magnetic field. The ring-shaped trapping potential 
can be created by a quadrupole storage ring trap [U EH , 
a linear rf multipole trap [Ijj, |20j or other methods. As 
an example, we consider N identical ions of mass M and 
charge q in a ring trap and a uniform magnetic field B 
(Fig. lb). The magnetic field is parallel to the axis of the 
trap. It is very weak so that it does not affect trapping. 
The equilibrium diameter of the ion ring is d. Figure 
lc shows examples of the trapping potentials for a 9 Be + 
ion in the radial plane of a quadrupole ring trap and a 
linear octupole trap, whose details are discussed in the 
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FIG. 1. Schematic of creating a space-time crystal, a, 

A possible structure of a space-time crystal. It has periodic 
structures in both space and time. The particles rotate in 
one direction even at the lowest energy state, b, Ultracold 
ions confined in a ring-shaped trapping potential in a weak 
magnetic field. The mass and charge of each ion are M and q, 
respectively. The diameter of the ion ring is d, and the mag- 
netic field is B. c, The pseudo-potentials (V ex t) for a 9 Be + ion 
in a quadrupole ring trap (solid curve) and a linear octupole 
trap (dashed curve) along the x or y axis. See Appendix for 
details. 



Appendix section. 

When the average kinetic energy of ions (ksT/2, where 
ks is the Boltzmann constant and T is the temperature) 
is much smaller than the typical Coulomb potential en- 
ergy between ions, i.e. T <C Nq 2 /(2n 2 eokBd), the ions 
form a Wigner ring crystal. For ions in a ring crystal, we 
can expand the Coulomb potential around equilibrium 
positions to the second order. So the many-body Hamil- 
tonian (see the Appendix section) becomes quadratic. 
We can diagonalize the quadratic Hamiltonian by intro- 
ducing a set of N normal coordinates qj and normal 
momenta p'y The normal coordinate and momentum 
of the collective rotation mode are q\ — -4= V\ 8 4 and 
■p' x = — L= pj , respectively. The remaining N — 1 nor- 
mal coordinates correspond to relative vibration modes. 
Choosing the potential energy at equilibrium positions 
as the origin of energy, the Hamiltonian of the system 



H = 
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N 
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(1) 

where h = h/(2ir) is the reduced Planck constant, 
a = qnd 2 B / (Ah) is the normalized magnetic flux, rj 2 — 
q 2 Md/{&irh 2 £Q), and uij is the normalized normal mode 
frequency. 

The eigenstate of this Hamiltonian is a product state 

tp[q) = e lk ^ J]^2 ^ife). where fjilj) U > 2 ) are ei S en - 
f unctions of harmonic oscillators. Resetting the ground 
state energy of the relative vibrations to be the origin of 
energy, the energy of this state is 



E = 



2h 2 
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where rij = 0, 1, 2, ■ ■ • (j > 2) are the occupation num- 
bers of the relative vibration modes. The wavefunction 
has to be symmetric with respect to the exchange of 
two identical bosonic ions, and has to be antisymmetric 
with respect to the exchange of two identical fermionic 
ions. If rij = for all j > 2, i.e. the relative vibration 
modes are in their ground states, the periodic bound- 
ary condition and symmetry property of the wavefunc- 
tion requires e lfe ( 27r /v / ^) = i for identical bosonic ions, 
and e l ' £ ( 27T 7v / A f ) _ (— i)^- 1 fo r identical fermionic ions. 
Thus for identical bosonic ions (e.g., 9 Be + ions), we have 
k = ni\/~N for all N, where n\ = 0, ±1, ±2, • • • . For iden- 
tical fermionic ions (e.g., 24 Mg + ions), we have k — nivN 
if A^ is an odd number and k = (n\ + \)yN if N is 
an even number. This results a qualitatively different 
rotation behavior for fermions and bosons that can be 
observed experimentally. 

The lowest normalized relative vibration frequency is 
u>2 = 2.48 when N = 10 and will increases as N increases 
2l. w 2 w v /0.32ATln(0.77AT) for large N. r) = 4.0 x 10 4 
for 9 Be + ions in a d = 10 /mi ring trap. So it costs a lot 
of energy to excite the relative vibration modes. Thus 
we have rij = for all j > 2 at lowest energy states. For 
an ion ring of identical bosonic ions, the energy E ni and 
the angular frequency u ni of the rii-th eigenstate of the 
collective rotation mode are: 

2/Vft 2 

^^(nr-a^^K-a) 2 , 



uj ni = uj*(ni - a) 



Md 2 
Ah 

(ni - a), 



(3) 



Md 2 



where E* = 2NH 2 /{Md 2 ) and u* = Ah/{Md 2 ) arc the 
characteristic energy and the characteristic frequency 
of the collective rotation, respectively. For identical 
fermionic ions, the results are the same as Eq. ([3]) if 
A^ is an odd number, and n\ should be changed to ri\ + 5 
if A^ is an even number. 
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FIG. 2. The energy levels and rotation frequencies of 
trapped ions, a, The energy levels of identical bosonic ions 
as a function of the magnetic flux a. The quantum number 
ni is labeled on each curve. The angular frequency of the 
persistent rotation as a function a is shown in b for an even 
number of fermionic ions, and c for bosonic ions, a and c are 
also applicable to an odd number of fermionic ions, d, The 
angular frequency of the persistent rotation of a bosonic ion 
ring in a constant magnetic field Bo > as a function of its 
normalized diameter d/do. 



In classical mechanics, the angular velocity of the low- 
est energy state is always lo/uj* = 0, which means that 
the ions do not rotate. In quantum mechanics, how- 
ever, u/uj* — is not an eigenvalue when the normalized 
magnetic flux a is not an integer or half of an integer. 
So the ions can rotate persistently even at the ground 
state. Since the ions are in the ground state already, 
there is no radiation loss due to the rotation. The rota- 
tion frequency is independent of the number of ions in the 
ring. The energy gap between the ground state and the 
first excited state is AE = Nh 2 /(Md 2 ) when a = 1/4. 
AE — > oo when N — > oo. Thus the persistent rotation 
of identical ions is a macroscopic quantum phenomenon 
and is robust for observation when N is large, which is 
important for a time crystal Q. This is very different 
from the situation of a rigid body with mass NM and 
charge Nq in a ring, for which the maximum energy gap 
between the ground state and the first excited state is 
AE rigld = 2h 2 /(NMd 2 ). AE ngid -> when N -> oo. 
When we consider a rigid body, we only minimize the en- 
ergy of its center-of-mass (c.o.m.) motion. Here we mini- 
mize the total energy of the whole system, which depends 
on the internal interaction and the exchange symmetry 
of the system. If the relative vibration modes are not 
in their ground states, the symmetry requirement of the 



c.o.m motion of identical particles is relaxed and the en- 
ergy gap between different c.o.m. motion states becomes 
smaller. The result becomes the same as that of a rigid 
body when all particles are different from each other. 

Figure 2a shows the lowest energy levels of an ion ring 
consisting identical bosonic ions. When the magnetic flux 
satisfies —0.5 < a < 0.5, the n\ — state is the ground 
state. As a increases above 0.5, n\ = 1 state becomes the 
ground state. Similar things happen whenever a acrosses 
half integer values. As a result, the angular frequency of 
the persistent rotation of the ground state is a periodic 
function of the magnetic flux (Fig. 2c). Figure 2d shows 
the rotation frequency of a bosnic ion ring in a constant 
positive magnetic field Bq as a function of its normalized 
diameter d/do, where do = y/ih/ (irqBo). The rotation 
frequency in Fig. 2d is normalized by Wg = qBo/2M. 
The ground state is m = when d/do < 1/V2- The 
rotation frequency is independent of the ring diameter 
and the number of ions in the ring for this state, and 
oscillates and decreases to when d/do increases above 
l/v2- If we confine many ions in a harmonic trap to 
form a 3D spatial crystal, ions in the crystal will rotate 
with the same angular frequency Wq and form a 4D space- 
time crystal when the outer diameter of the ion crystal is 
smaller than do /V%- If we confine ions in two concentric 
ring traps with diameters larger than do/V%, the rotation 
frequencies of the two rings can be different or the same, 
depending on the interaction between ions in different 
rings. When the ratio of the rotation frequencies of the 
two rings is an irrational number, the ions have an order 
in time but cannot reproduce their positions simultane- 
ously in a finite period. Thus we have a time quasicrystal, 
in analog of a conventional spatial quasicrystal 22 1. 

The persistent rotation of trapped ions can be detected 
by measuring the current generated by the ions, or in- 
ferred by probing the energy levels of the ion ring. More 
importantly, we can observe the persistent rotation di- 
rectly by measuring the ion positions twice when N is 
large. For example, if we have an ion ring consisting N 
identical 9 Be + ions, we can first use a pulse of two co- 
propagating laser beams to change the hyperfine state of 
one (or a small fraction) of the ions by stimulated Ra- 
man transition and use this ion as a mark (qubit memory 
coherence time greater than 10 s has been demonstrated 
with 9 Be + ions [23j). Both laser beams are parallel to the 
axis of the ion ring and have waists of wq- We assume 
that wo is larger than the separation between neighbor- 
ing ions and the pulse is very weak so that on average 
less than one ion is marked. This two-photon process lo- 
calized the position of the mark ion with an uncertainty 
of about Ax ~ The recoil momentum of the ion 

due to a stimulated Raman transition is p\ — P2, where 
pi and p2 are the momenta of photons in the two laser 
beams. The amplitude of the transverse momentum of 
each photon in a Gaussian beam with waist of wo is about 
h/uio- Thus the momentum of the ion ring is changed by 



4 



13 



b 10" 1 




gular frequency of the ions is 



"ICT-i 



T 2.6x1 6 




2.6x10° 



10 
d (urn) 



FIG. 3. The temperature dependence of the persis- 
tent rotation of trapped ions, a, The average angular 
frequency of the persistent rotation of identical bosonic ions 
as a function of the temperature. From top down, the mag- 
netic flux increases from -0.45 to 0.45. b, The characteris- 
tic temperature (left axis) and the characteristic frequency 
(right axis) of the persistent rotation of an ion (or electron) 
ring consisting fOO identical ions (or electrons) as a function 
of the diameter. From top down, the trapped particles are 
electrons, 9 Be + ions, and 24 Mg + ions. 



about Ap ps y/^h/iVQ. Ap should be smaller than the 
absolute value of the initial momentum of the ion ring, 
which is Nh/(2d) when a — 1/4. Thus the waists of 
lasers need to satisfy 2y2ci/iV < wo < V%d in order 
to localize the position of an ion without significantly 
alter the initial momentum of the ion ring. This condi- 
tion can be fulfilled when N is large. Then we can use a 
global probe laser which is only scattered by the mark ion 
[24J (state-dependent fluorescence) to measure its angu- 
lar displacement (A9 m k) after a time separation At. The 
displacement of the mark is about A9 m k ps w*(ni — a) At 
when TV is large. The mark ion repeats its position when 
At fa 2irl/[uj*(ni ~ a)] , where I is an integer. After the 
measurement, we can cool the ions back to the ground 
state and repeat the experiment again. 
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n i — — oo 



-E n Jk B T 
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where the partition function is Z — Y^m e~ B "i/ feBT '. ft 
is convenient to define T* = E*/k B — Nhio*/2k B as 
the characteristic temperature for the ring of ions. T* 
increases when N increases. 

Figure 3a shows the average rotation frequency of an 
ion ring consisting identical bosonic ions as a function of 
the temperature. At very low temperatures (T <C T*), 
the average rotation frequency is independent of the tem- 
perature. As the temperature increases, the probability 
of the ion ring occupying excited states increases. These 
states have positive and negative velocities alternately, 
which cancel each other. As a result, the amplitude of the 
average rotation frequency drops to zero at high temper- 
atures (T > T*). So T* can be considered as the phase 
transition temperature of the space-time crystal. T* and 
uj* as a function of the diameter of an ion ring (or an elec- 
tron ring) are displayed in Fig. 3b. For a d = 100 /im ion 
ring consisting 100 9 Be+ ions, T* = 1.1 nK and u* = 2.8 
rad/s. T* is larger for smaller ion rings. 

In order to experimentally realize such a space-time 
crystal with trapped ions, we need to confine ions tightly 
to have a small d, and cool the ions to a very low temper- 
ature. Recently, cylindrical ion traps with inner radius 
as small as 1 /im have been fabricated [25j. Simulations 
suggested that it is also possible to confine charged par- 
ticles with a nanoscale rf trap [26| . The challenge is that 
ions must be cooled to below 1 /^K for a microscale trap 
(Fig. 3b). We propose to first add a pinning poten- 
tial to confine ions with MHz trapping frequencies in the 
circumference direction. A combination of Doppler cool- 
ing and resolved-sideband cooling can be used to cool 
the ions to the ground state of the MHz trap [l5|, 27 1. 
The system is in the ground state of the ring-shaped 
trapping potential after ramping down the pinning po- 
tential adiabatically. For T* = 1.1 nK, the ramping 
down time should be longer than 7 ms. Another way 
to achieve an ultralow temperature of ions is to put the 
ions near a Bosc-Einstcin condensate of neutral atoms 
[28j . which has been cooled to below 0.5 nK by adiabatic 
decompression [2^| . 

In conclusion, we propose a method to create and ob- 
serve a space-time crystal experimentally with trapped 
ions. We also discuss about how to create a time qua- 
sicrystal. The space-time crystals of trapped ions provide 
a new dimension for studying many-body physics and 
may have potential applications in quantum information 
for simulating other novel states of matter [301 ] . 



To study the effects of finite temperatures on persis- 
tent rotation, we assume the trapped ions are at thermal 
equilibrium with temperature T. Then the average an- 
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APPENDIX 

Quadrupole ring trap. A rf quadrupole ring trap 
can be formed by many segments of quadrupole elec- 
trodes 17 1 or four ring electrodes [18(. The harmonic 
pseudo-potential for a single ion is V ex t(r, z) = ^Muj 2 (r— 

r m in) 2 + \Muj 2 z 2 , where r = \/ x 2 + y 2 is the radial po- 



where pj 



;2fi a 
d d9j 



is the canonical momentum opera- 



sition of the ion, r m in is the radius of the ring trap, u> r 
and w z are trapping frequencies in the radial and axial 
(z) directions, respectively. If the inner radius of the trap 
(ro , which is half of the separation between opposite elec- 
trodes [17]) is much smaller than r m i n , the trapping fre- 
quencies are approximately ui r = ui z = qV r f/(^/2Milr 2 ). 
Here V r f is the amplitude of the rf voltage applied to 
adjacent electrodes, ri/27r is frequency of the rf voltage. 
In Fig. lc, uj r — 2tt x 5 MHz and r m i n = 20 /im for the 
quadrupole ring trap. 

Linear rf multiple trap. The pseudo-potential for 
a single ion in a linear rf multipole trap is cylindrically 
" 2C 



symmetric 



1!: 



V ext (r,z) 



k 2 q 2 V 2 



rf 



(-) 



2k-2 



PqVdc , 2 _ r \ 
.2 \ Z o )> 



where 2k is the number of poles of the trap, /3 is the 
geometrical factor of the trap, Vdc is the static voltage 
applied to the outside sets of segmented electrodes, and 
2zo is the length of the center part of the segmented 
trap. The trapping frequency in the axial direction is 
uj z = \J 2(3qVdc/{MzQ) . The trapping potential has its 
minimum value in the radial direction at 



= r 



2MOcj 



z> 



qV r fk\/k — 1 



l/(fc-2) 



The trapping frequency in the radial direction is u r = 
Vfc — 2w 2 near r m i n . Fig. lc shows an example of the 
effective trapping potential for a 9 Be + ion in the radial 
plane of an octupole trap (2k = 8). The parameters used 
for the calculation of Fig. lc are ro = 200 /im, V r f = 200 
V, w 2 /2tt = 0.5 MHz, and fL/2ir = 94 MHz. 

Many-body Hamiltonian. As we are interested in 
the rotation of ions in a ring at very low temperatures, 
our Hamiltonian only involves the angular coordinate 6j 
(j = 1, • • • , N) of trapped ions: 



H 



N 

£ 



(Pj - gA e f 
2M 



V q - 

W|sin(^)| 



tor, Ag = Bd/A = /p/ijrd) is the vector potential, and eo 
is the electric constant of vacuum. cf> = ird 2 B/4 is the 
magnetic flux inside the ion ring. 
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